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Abstract. This paper describes an efficient and unique 
method for computing the shortwave direct radiative ef- 
fect (DRE) of aerosol residing above low-level liquid-phase 
clouds using CALIOP (Cloud-Aerosol Lidar with Orthogo- 
nal Polarization) and MODIS (Moderate Resolution Imaging 
Spectroradiometer) data. It addresses the overlap of aerosol 
and cloud rigorously by utilizing the joint histogram of cloud 
optical depth and cloud top pressure while also accounting 
for subgrid-scale variations of aerosols. The method is com- 
putationally efficient because of its use of grid-level cloud 
and aerosol statistics, instead of pixel-level products, and a 
precomputed look-up table based on radiative transfer cal- 
culations. We verify that for smoke and polluted dust over 
the southeastern Atlantic Ocean the method yields a sea- 
sonal mean instantaneous (approximately 13:30 local time) 
shortwave DRE of above-cloud aerosol (ACA) that gener- 
ally agrees with a more rigorous pixel-level computation 
within 4 %. We also estimate the impact of potential CALIOP 
aerosol optical depth (AOD) retrieval bias of ACA on DRE. 
We find that the regional and seasonal mean instantaneous 
DRE of ACA over southeastern Atlantic Ocean would in- 
crease, from the original value of 6.4 W m -2 based on oper- 
ational CALIOP AOD to 9.6 Wm -2 if CALIOP AOD re- 
trievals are biased low by a factor of 1.5 (Meyer et al., 
2013) and further to 30.9 W m -2 if CALIOP AOD retrievals 
are biased low by a factor of 5 as suggested in Jethva 
et al. (2014). In contrast, the instantaneous ACA radiative 


forcing efficiency (RFE) remains relatively invariant in all 
cases at about 53 Wm -2 AOD -1 , suggesting a near-linear 
relation between the instantaneous RFE and AOD. We also 
compute the annual mean instantaneous shortwave DRE of 
light-absorbing aerosols (i.e., smoke and polluted dust) over 
global oceans based on 4 years of CALIOP and MODIS 
data. We find that given an above-cloud aerosol type the op- 
tical depth of the underlying clouds plays a larger role than 
above-cloud AOD in the variability of the annual mean short- 
wave DRE of above-cloud light-absorbing aerosol. While we 
demonstrate our method using CALIOP and MODIS data, it 
can also be extended to other satellite data sets. 


1 Introduction 

The shortwave direct radiative effect (DRE) of aerosols at 
the top of the atmosphere (TOA) is strongly dependent on 
the reflectance of the underlying surface. Over dark surfaces 

(e.g., ocean, vegetated land), the scattering effect of aerosols 

is generally dominant, leading to negative DRE (i.e., cooling) 
at TOA (Yu et al., 2006). In contrast, when light-absorbing 
aerosols occur above clouds or other bright surfaces (such as 
snow, ice, and desert), aerosol absorption is significantly am- 
plified by cloud or surface reflection, offsetting or even ex- 
ceeding the scattering effect of the aerosol, leading to a less 
negative or even positive (i.e., warming) TOA DRE (Abel et 
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al., 2005; Keil and Haywood, 2003; Twomey, 1977). There- 
fore, in order to understand the full complexity of aerosol 
radiative effects on climate, it is important to quantify the 
DRE under both clear-sky and cloudy-sky conditions. The 
DRE of aerosols in clear-sky regions has been extensively 
studied and is relatively well constrained based on advanced 
satellite remote sensing measurements acquired in the last 
decade (Yu et al., 2006). However, current model simulations 
show a large intermodel spread in cloudy-sky DRE (Schulz et 
al., 2006), which results from intermodel differences in both 
aerosol and cloud properties (Schulz et al., 2006; Stier et al., 
2013). Therefore, there is a clear need for an observational 
constraint on the DRE of above-cloud aerosol (ACA). 

Recent advances in satellite remote sensing techniques 
have provided an unprecedented opportunity for studying 
the DRE of ACA. In particular, the availability of mea- 
surements from the space-borne Cloud-Aerosol Lidar with 
Orthogonal Polarization (CALIOP) sensor onboard NASA’s 
Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Ob- 
servations (CALIPSO) satellite has provided a revolution- 
ary global view of the vertical distribution of aerosols and 
clouds (e.g., Winker et al., 2013). Using CALIOP aerosol 
and cloud layer products, Devasthale and Thomas (2011) 
found frequent occurrences of aerosols residing above low- 
level clouds in several regions of the globe. In particular, they 
found a high frequency of smoke occurrence over low clouds 
in the southeastern Atlantic, western coasts of South Amer- 
ica (e.g., Colombia, Ecuador, and Pern) and southern Asia. 
These authors also found a high frequency of natural and 
polluted dust aerosols overlapping low clouds off the western 
coast of Saharan Africa in boreal summer and over boundary 
layer clouds in the eastern coast of China in boreal spring 
(see Fig. 3 of Devasthale and Thomas, 2011). 

CALIOP measurements of ACA properties, in combina- 
tion with satellite cloud products from, for example, the 
Moderate Resolution Imaging Spectroradiometer (MODIS), 
have been used in several recent studies to derive the DRE 
of ACA with radiative transfer simulations (e.g., Chand et 
al., 2009; Costantino and Breon, 2013b; Meyer et al., 2013; 
Oikawa et al., 2013). Chand et al. (2009) used CALIOP 
above-cloud AOD retrievals (Chand et al., 2008) and Terra- 
MOD1S cloud products, both aggregated to 5° gridded 
monthly means, to calculate the radiative effects of smoke 
transported above the low-level stratocumulus deck in the 
southeastern Atlantic. The spatiotemporal aggregation of 
both CALIOP and MODIS data to coarse-gridded monthly 
means obscures the potential influence of cloud and aerosol 
variability on the DRE. In particular, using the grid-box mean 
cloud optical depth for DRE calculation might lead to biases 
in DRE due to the plane-parallel albedo bias (Oreopoulos et 
al., 2007). Moreover, the MODIS level-3 aggregation algo- 
rithm samples all liquid water clouds, regardless of whether 
or not there is an aerosol layer above clouds. As a result, 
the total population of liquid water clouds in the MODIS 
level-3 products (daily or monthly) may be significantly 


Z. Zhang et al.: Computing DRE of above cloud aerosols 

different from that of the below-aerosol-only cloud popula- 
tion. Therefore, using level-3 MODIS products without dis- 
tinguishing below-aerosol-only from total cloud population 
can potentially lead to significant errors. The problem could 
be further complicated by biases in MODIS cloud retrievals 
associated with the presence of overlying light-absorbing 
aerosols. When a cloud pixel is contaminated by overly- 
ing light-absorbing aerosols the MODIS cloud optical depth 
(COD) retrieval is generally biased low (e.g., Coddington et 
al., 2010; Haywood et al., 2004; Jethva et al., 2013; Wilcox, 
2010), an effect not considered in most previous studies (e.g., 
Chand et al., 2009; Costantino and Breon, 2013b; Oikawa et 
al., 2013). Most recently however, Meyer et al. (2013) collo- 
cated CALIOP above-cloud AOD and Aqua-MODIS cloud 
properties at the pixel level, and the DRE was then computed 
at these individual collocated pixels. They found that correct- 
ing the MODIS COD bias due to ACA contamination leads 
to a more positive ACA DRE. Such rigorous collocation has 
obvious advantages as it takes into account the subgrid vari- 
ability of clouds and aerosols, but is however computation- 
ally expensive since it requires large amounts of pixel-level 
data that make global scale and multiyear studies challeng- 
ing. 

The objective of this paper is to describe a novel method 
for computing the DRE of ACA. This method attempts to 
balance the need for computational efficiency with the need 
for rigorous treatment of aerosol-cloud overlap and small- 
scale variability of aerosol and clouds. Our method has sev- 
eral unique features: (1) it takes subgrid-scale cloud and 
aerosol variation into account in DRE computations, (2) it 
treats the overlap of aerosol and cloud rigorously by utilizing 
the joint histogram of COD and cloud top pressure (CTP) in 
the MODIS level-3 product, and (3) it is computationally ef- 
ficient because of the use of a precomputed look-up table of 
ACA DRE. 

In the following sections, we briefly introduce the 
CALIOP and MODIS data used (Sect. 2), describe the key 
assumptions and features of the novel method (Sect. 3), val- 
idate it through comparison with pixel-level computations as 
in Meyer et al. (2013) (Sect. 4), and conclude with a sum- 
mary and discussion (Sect. 5). 


2 Satellite data 

In Meyer et al. (2013), the MODIS level-2 cloud product 
is collocated with the CALIOP level-2 aerosol product for 
every pixel along the CALIOP track and the computation 
of instantaneous DRE is performed pixel by pixel. Then, 
the pixel-level DRE results are aggregated on a latitude- 
longitude grid for climatological study. If only the grid- 
level DRE is of interest, the pixel-by-pixel computation of 
DRE may not be efficient because of redundant computa- 
tions. For example, if two pixels with the same above -cloud 
AOD and below-cloud COD occur within the same grid box, 
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they evidently have the same ACA DRE, but the radiative 
transfer computation is nevertheless performed twice in the 
pixel-by-pixel method. As shown in Sect. 3, a more efficient 
way is to compute the DRE statistically using the probabil- 
ity density function (PDF) of above-cloud AOD and below- 
cloud COD. In this study, we use the CALIOP level-2 aerosol 
and cloud layer product (V3.01) to derive the statistics of 
ACA properties and the MOD1S level-3 daily cloud prod- 
uct for cloud property statistics. It is important to note that 
our method is not limited to CALIOP and MOD1S prod- 
ucts, but is also applicable to other satellite data sets, such as 
above-cloud aerosol retrievals from POLDER (POLarization 
and Directionality of the Earth’s Reflectances) (Waquet et 
al., 2009) and OMI (Ozone Monitoring Instrument) (Torres 
et al., 2012), and cloud retrievals from ISCCP (International 
Satellite Cloud Climatology Project) (Rossow and Schiffer, 
1999) and SEVIRI (Spinning Enhanced Visible and Infrared 
Imager) (Schulz et al., 2009). 

2.1 CALIOP level-2 aerosol and cloud layer products 

Since its launch in 2006, the space-borne lidar CALIOP 
has continuously acquired, with near global (albeit instan- 
taneously sparse) coverage, attenuated backscatter measure- 
ments at 532 and 1064 nm, including linear depolarization 
information at 532 mn (Winker et al., 2009). The CALIOP 
level-2 retrieval algorithm consists of several steps. First, 
a “feature finder” algorithm and cloud-aerosol discrimina- 
tion (CAD) algorithm are used to detect aerosol and cloud 
layers, and record their top and bottom heights and layer- 
integrated properties (Vaughan et al., 2009). Second, the de- 
tected aerosol layers are further classified into six subtypes 
(i.e., polluted continental, biomass burning, desert dust, pol- 
luted dust, clean continental and marine) (Omar et al., 2009) 
and the detected cloud layers are assigned different thermo- 
dynamic phases (Hu et al., 2007a) based on the observed 
backscatter, color ratio and depolarization ratio. Third, a pri- 
ori lidar ratios, preselected based on aerosol subtype and 
cloud phase, are used to derive the extinction of an aerosol 
or cloud layer from the attenuated backscatter profile (Young 
and Vaughan, 2008). 

In this study, we use CALIOP level-2 version 3.01 
aerosol and cloud layer products at a nominal 5 km 
horizontal resolution (i.e., CAL_LID_L2_05kmALay and 
CAL_LID_L2_05kmCLay) for aerosol-cloud overlap detec- 
tion, and for information on aerosol layer properties, includ- 
ing type, aerosol optical depth (AOD), and layer top and bot- 
tom height. In addition to physical properties, the CALIOP 
layer products also provide various metrics and flags on 
data quality assurance. These include CAD score (Liu et al., 
2009), horizontal averaging scale, extinction quality control 
(QC) flag, and estimated uncertainty of layer AOD. In this 
study, we apply these metrics following best practices pro- 
vided by the CALIPSO science team to screen for reliable 
retrievals (e.g., Winker et al., 2013) (see Table 1). 
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Table 1 . Quality control metrics used for screening the CALIOP 
aerosol layer product. 


Criterion 

CAD score < — 30 

Horizontal averaging < 20 km 

Extinction_QC_532 0 or 1 

Feature_optical_depth_uncertainty_532 < —99.5 


It should be noted here that the current version of CALIOP 
operational aerosol retrieval algorithm (V3.01) appears to 
significantly underestimate the AOD of above-cloud aerosol 
layer according to recent studies (Jethva et al., 2014; Kacene- 
lenbogen et al., 2013, 2011; Liu et al., 2013; Torres et al„ 
2013; Waquet et al., 2013b). The main reason is that after 
strong attenuation by the upper part of an aerosol layer, the 
532 nm attenuated backscatter of the lower part of aerosol 
layer is often too small (Torres et al., 2013). As a result, the 
current CALIOP feature detection algorithm often cannot re- 
solve the full depth of the aerosol layer, leading to low biases 
in retrieved AOD (Jethva et al., 2014; Liu et al., 2013; Torres 
et al., 2013). This issue is more acute for smoke than dust 
because smoke has a much larger lidar ratio. At the moment, 
the CALIOP operational team is investigating the possibility 
of using the algorithm described in Chand et al. (2008) and 
Hu et al. (2007b) for ACA retrievals (Liu et al., 2013). This 
alternate method utilizes the reflected lidar signal from the 
bright cloud layer underneath to derive the two-way trans- 
mittance and thereby the AOD of the ACA layer. Because 
the backscatter of a cloud layer is usually very strong, the 
two-way transmittance method is less affected by the strong 
attenuation of the ACA layer and is therefore expected to 
alleviate the aforementioned problem. Lidar based AOD re- 
trievals are also known to suffer from other issues, such as 
the background solar noise during daytime. These issues are 
beyond the scope of this study, but are nevertheless discussed 
in the uncertainty analysis of Sect. 3.4. 

In addition to retrieval errors and uncertainties, another 
limitation of CALIOP data is the small sampling rate (i.e., 
only along track). In order to compute the DRE of ACA 
over a given latitude-longitude grid box, we assume that 
the aerosol property statistics retrieved by CALIOP along 
its narrow track represent the statistics over the whole grid 
box, i.e., that AOD PDFs are identical. This assumption con- 
stitutes an uncertainty in our DRE computation. Due to a 
lack of satellite -based wide-swath ACA data sets, however, 
it is difficult to determine the size of this uncertainty. Re- 
cently, several novel methods have been developed to retrieve 
ACA properties from passive sensor observations (Jethva et 
al., 2013; Torres et al., 2012; Waquet et al., 2009, 2013a), 
which will help improve our understanding of the subgrid 
ACA variability when they become available to public. 
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Finally, we emphasize two more points. First, none of the 
aforementioned problems with CALIOP data, e.g., smoke 
AOD bias, retrieval uncertainties, and small sampling rate, 
are unique to our method. Any method that uses CALIOP 
data faces the same challenges. Second, our method is not 
limited only to CALIOP data. We choose to use the CALIOP 
product in this study solely because it is the only publicly 
available ACA product at the moment. Our method can also 
be applied to other ACA retrieval products based on, for ex- 
ample, POLDER (Waquet et al., 2009), MODIS (Jethva et 
ah, 2013), and OMI (Torres et ah, 2012) observations when 
they become available to the public. In fact, as discussed 
later, the advantage of our method in terms of computational 
efficiency is even greater when applied on retrievals from 
passive sensors. 

2.2 MODIS daily level-3 cloud property product 

This study computes the grid-level ACA DRE using the 
statistics of aerosol and cloud properties, instead of pixel- 
by-pixel computation as in Meyer et ah (2013). We use the 
collection 5 (C5) Aqua-MODIS level-3 daily gridded atmo- 
sphere product MYD08 D3 for the statistics of cloud proper- 
ties and other parameters, such as solar zenith angle, needed 
for ACA DRE computations. 

The MODIS level-3 (i.e., grid level) product contains 
statistics computed from a set of level-2 (i.e., pixel level) 
MODIS granules. As summarized in Platnick et ah (2003), 
the operational level-2 MODIS cloud product provides cloud 
masking (Ackerman et ah, 1998), cloud top height retrieval 
based on CO 2 slicing or the infrared window method (Men- 
zel et ah, 1983), cloud top thermodynamic phase determina- 
tion (Menzel et ah, 2006), and cloud optical and microphysi- 
cal property retrieval based on the bispectral solar reflectance 
method (Nakajima and King, 1990). In addition to these 
cloud parameters, the level-2 products also provide pixel- 
level runtime quality assessment (QA) information, which 
includes product quality as well as processing path informa- 
tion. All MODIS level-2 atmosphere products, including the 
cloud, aerosol and water vapor products, are aggregated to 
1° spatial resolution on a daily (product name MYD08 D3 
for Aqua MODIS), eight-day (MYD08 E3), and monthly 
(MYD08 E3) basis. Aggregations include a variety of scalar 
statistical information (mean, standard deviation, max/min 
occurrences) and histograms (marginal and joint). A partic- 
ularly useful level-3 cloud product for this study is the daily 
joint histogram of COD vs. CTP, derived using daily counts 
of successful daytime level-2 pixel retrievals that fall into 
each joint COD-CTP bin. Eleven COD bins, ranging from 
0 to 100, and 13 CTP bins, ranging from 200 to 1000 mb, 
comprise the histogram. As discussed below, the COD-CTP 
joint histogram allows for identification of the portion of the 
cloud population that lies beneath the aerosol layer found 
by CALIOP, as well as the corresponding COD probability 
distribution needed for DRE estimation. In addition to the 
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COD-CTP joint histogram, we also use the gridded mean 
solar and sensor zenith angles for calculating DRE and cor- 
recting the COD bias due to the presence of ACA. 

It should be noted that the level-3 daily product 
MYD08 D3 contains statistics computed from a set of level- 
2 MODIS granules that theoretically span a 24 h interval 
(Hubanks et al., 2008). However, for cloud parameters re- 
trieved only during daytime, such as COD and cloud droplet 
effective radius (CER), only daytime level-2 files are used to 
compute the level-3 daily statistics. These are called daytime 
only SDSs (scientific data sets) in level-3 products. Strictly 
speaking, the daytime only SDSs of only those 1° grid cells 
between approximately 23° N and 23° S come from a sin- 
gle MODIS overpass. The tropical southeastern Atlantic re- 
gion, where transported smoke aerosols are often observed 
above low-level stratocumulus clouds, is well within this 
range (about 10°N-30°S see Fig. 3). The COD statistics 
in the MYD08 D3 product for this region are therefore de- 
rived from a single Aqua-MODIS overpass that can be col- 
located with CALIOP observations (see Sect. 3.1 for de- 
tails on collocation). The DRE computed based on the collo- 
cated data set is therefore instantaneous DRE at Aqua cross- 
ing time (13:30 LT, local time) that are directly comparable 
to the pixel-by-pixel results in Meyer et al. (2013). Pole- 
ward of 23°, MYD08 D3 statistics are derived from averag- 
ing several overlapping orbits approximately 100 min apart 
(Hubanks et al., 2008). As a result, strictly speaking the 
DRE computed for mid- and high-latitude regions based on 
MYD08 D3 data is not instantaneous DRE. We emphasize 
that this is not a limitation of our method, but an inher- 
ent characteristic of the MODIS level-3 product. Finally, it 
needs to be noted that MODIS cloud property products are 
subject to retrieval uncertainties caused by various factors 
such as subpixel inhomogeneity, 3-D radiative transfer ef- 
fects and presence of precipitation in clouds (Zhang, 2013; 
Zhang and Platnick, 2011; Zhang et al., 2012). The uncer- 
tainties in cloud property retrieval could translate to uncer- 
tainty in the ACA DRE computation, which will be investi- 
gated in future studies. 


3 Methodology 
3.1 Theoretical basis 

As in previous investigations (e.g., Chand et al., 2008, 2009; 
Costantino and Breon, 2013b; Meyer et al., 2013), we fo- 
cus on the simplest case of overlapping aerosol and cloud, 
i.e., a single layer of aerosol overlying a single layer of low- 
level liquid-phase clouds, which is commonly observed in 
many regions of the globe (Devasthale and Thomas, 2011). 
More complex situations certainly exist, such as an aerosol 
layer located in between high and low clouds, or an aerosol 
layer overlying multiple layers of clouds. However, identi- 
fication of such situations are either beyond the detection 
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Figure 1 . A schematic example to illustrate how CALIOP aerosol 
layer height information is used in our method to determine the pop- 
ulation of liquid-phase clouds below the aerosol layer in the MODIS 
COD-CTP joint histogram. 


capabilities of CALIOP or relatively rare (Devasthale and 
Thomas, 2011). As such, they are not considered here and 
left for future research. 

To illustrate the theoretical foundation of the method, con- 
sider the schematic example in Fig. 1. For a given grid box 
(e.g., 1° x 1° in case of MODIS level-3 data), the gridded 
mean instantaneous broadband shortwave DRE ( (DRE) A(:A ) 
averaged over all ACA pixels within the grid box is given as 

OO OO 

(DRE) aca = J J DRE(r c ,r a )p(r c ,r a )t/r c t/r a , (1) 

o o 

where p( r c , r a ) is the joint PDF of the above-cloud AOD at 
532 nm (r a ) and below-aerosol COD (r c ) of ACA pixels. We 
note that, in addition to r a , DRE also depends on the spectral 
variation of aerosol and cloud optical depth, spectral single 
scattering albedo and asymmetry factor, wavelength depen- 
dencies not explicitly shown in this equation. These proper- 
ties are computed using a Mie scattering code (Wiscombe, 
1980) based on the aerosol model described in Meyer et 
al. (2013). The dependencies on solar zenith angle, surface 
reflectance, cloud particle effective radius, and atmospheric 
profile are also omitted from the equation; solar zenith angle 
and surface reflectance are expected to have only minor vari- 
ation within the grid box, while the impact of cloud particle 
effective radius and atmospheric profile on shortwave DRE is 
relatively small. Since p(r c , r a ) describes the covariation of 
aerosols and clouds for the ACA pixels, it should ideally be 
derived from collocated CALIOP aerosol and MODIS cloud 
retrievals at pixel level as in Meyer et al. (2013). However, 
this requires large amounts of pixel-level data. Therefore, 
pixel-level collocation and radiative transfer simulation are 
too computationally expensive and cumbersome for multi- 
year global studies. 

A key assumption in our method, which allows us to 
avoid tedious pixel-level collocation, is that the subgrid 
level instantaneous spatial distribution of above-cloud AOD 
is statistically independent from the subgrid level instanta- 
neous spatial distribution of below-aerosol COD. Under this 


assumption, p( r c , r a ) = p( r c ) • p(r a ) and Eq. (1) reduces to 


(DRE) aca — 


// DRE(t c , r a )p(r c )dx c 
o Lo 


p(x a )dx a , (2) 


where p( r c ) and p( r a ) are the PDF of instantaneous below- 
aerosol COD r c and above-cloud AODr a , respectively, of 
ACA pixels. The advantage of Eq. (2) is that it allows p( r c ) 
and p ( r a ) to be derived separately and independently. This 
assumption is reasonable considering that transported ACAs 
and low-level boundary layer clouds are usually well sep- 
arated vertically (Devasthale and Thomas, 2011) and con- 
trolled by different meteorological conditions. The poten- 
tial coupling between the two is that overlying absorbing 
aerosols could influence the evolution of clouds through 
changing atmospheric stratification (Wilcox, 2010). How- 
ever, a recent observational study (Costantino and Breon, 
2013a) found no correlation between above-cloud AOD and 
below-aerosol COD, although correlations are found be- 
tween AOD and cloud droplet effect radius, as well as liq- 
uid water path. Moreover, it is important to stress that our 
assumption is that the instantaneous above-cloud AOD and 
below-aerosol COD are independent at subgrid scale. This 
assumption does not rule out the possibility that AOD and 
COD could be correlated at longer temporal (e.g., seasonal) 
and/or larger spatial (e.g., regional) scales through the ther- 
modynamic and radiative coupling (Wilcox, 2010, 2012). Fi- 
nally, as shown in Sect. 4, when we compare the DRE derived 
from pixel-level collocation (i.e., based on Eq. 1) with that 
from independent sampling of p{r c ) and p(x a ) (i.e., based 
on Eq. 2) the agreement is very good. 

In our method, the PDF of above-cloud AOD p( r a ) is de- 
rived from the CALIOP 5 km aerosol and cloud layer prod- 
ucts through the following steps: (1) for each 5 km CALIOP 
profile that falls within a given latitude-longitude grid box, 
we first search for an aerosol layer; (2) if an aerosol layer is 
detected and the quality metrics pass the quality assurance 
criteria summarized in Table 1, we then proceed to check 
for the presence of an underlying liquid-phase cloud layer 
within the profile using the CALIOP cloud layer product; 
and (3) if a cloud layer is present, the AOD of the aerosol 
layer is recorded for the derivation of the p( r a ) of the grid 
box. The bottom height of the aerosol layer is also recorded 
to derive the grid mean aerosol layer bottom height. Once all 
of the CALIOP profiles within the grid box are processed, we 
obtain the PDF of the above-cloud AODp(r a ) and the mean 
aerosol layer bottom pressure (Pbottom)- 

As schematically illustrated in Fig. 1, the PDF of below- 
aerosol COD p(x c ) is derived from the joint histogram of 
cloud optical depth and cloud top pressure (COD-CTP joint 
histogram) in the MODIS daily level-3 product, using the 
grid-mean aerosol-layer bottom pressure (Pbottom) derived 
above. For a given grid box, we first identify the popula- 
tion of liquid-phase clouds below the pressure level (^bottom}- 
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This subset, together with the AOD PDF/?(r a ), is then used 
to calculate DRE according to Eq. (2). 

In this study, we focus on the computation of instanta- 
neous DRE. To obtain diurnally averaged DRE, technically 
speaking one would simply need to integrate over time the 
instantaneous DRE. However, it is important to note that, in 
addition to diurnal variation of solar zenith angle, aerosol 
and cloud properties may also have significant diurnal cy- 
cles. In fact, it is known that the low cloud fraction over stra- 
tocumulus regimes, such as the southeastern Atlantic region, 
have a strong diurnal cycle (15-35 % of diurnal mean value) 
driven by cloud solar absorption (Wood et ah, 2002). A re- 
cent study by Min and Zhang (2014) indicates that using a 
constant cloud fraction based on Aqua-MODIS observations 
tends to result in significantly underestimated diurnal mean 
DRE even if the diurnal variation of solar zenith angle is 
considered in the computation. Therefore, the diurnal vari- 
ation of cloud properties is an important factor to consider in 
the intercomparison of DRE computations based on different 
data sets and intercomparison between observational study 
and modeling results. 

3.2 DRE look-up tables 

To speed up calculations, we use precomputed aerosol-type 
specific look-up-tables (LUTs), instead of online radiative 
transfer computation, when deriving the DRE of ACA. The 
concept of our LUTs is somewhat similar to the “radiative 
kernels” described in Hartmann et al. (2001) and Zelinka 
et ah (2012) for computing cloud radiative feedbacks. The 
LUT for each aerosol type consists of DREs at both TOA 
and surface (not used in this study) for various combinations 
of AOD, COD, CTP and solar zenith conditions. As such, 
once the aerosol type and AOD are known from CALIOP and 
COD, CTP and solar zenith angle are known from MODIS, 
and the corresponding DRE can be obtained through LUT 
interpolation. Note that the CALIOP only provides AOD at 
lidar wavelengths (e.g., 532 and 1064nm) for each aerosol 
type. Therefore, radiative transfer model-appropriate nar- 
rowband aerosol scattering properties, namely AOD, single- 
scattering albedo and asymmetry factor, are needed for the 
development of the LUT. The current version of LUT focuses 
on light-absorbing aerosols (e.g., smoke and polluted dust). 
In order to validate our method with more rigorous pixel- 
level computations, we adopt the narrowband aerosol optical 
properties of Meyer et ah (2013), who used the same radia- 
tive transfer code, in the computation of the current LUT. 
The aerosol model in Meyer et ah (2013) is extended from 
an absorbing aerosol model developed for the MODIS Col- 
lection 5 aerosol product (MOD04) (see Table 4 of Levy 
et ah, 2009). The MOD04 aerosol models define aerosol 
size distributions and refractive indices based solely on pre- 
scribed AOD at 550 mil (MODIS band 4; note that the ab- 
sorbing aerosol model used here assumes a constant index 
of refraction, 1.51-0.02/, at all wavelengths). At AOD = 0.5 
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Figure 2. A schematic illustration of our fast scheme to correct the 
COD retrieval bias in the MODIS cloud product due to overlying 
aerosol contamination. 


(550 mn), the single-scattering albedo of this model is about 
0.9 over the visible spectral region (see Fig. 7 of Meyer et 
ah, 2013), which is in the range of previously reported val- 
ues (e.g., Keil and Haywood, 2003; Myhre et ah, 2003). 
The current AOD bins (at 550 nm) in the LUT range from 
0.05 to 1.5, which covers most of the above-cloud AOD ob- 
served by CALIOP. The current COD bins, logarithmically 
spaced, range from 0.1 to 300. Following the MODIS level- 
3 data, the thirteen CTP bins range from 1000 to 200 mb. 
The solar zenith angle bins range from 0 to 80°. Radiative 
transfer computations are carried out using the RRTM-SW 
(rapid radiative transfer model, shortwave) model (Clough et 
ah, 2005; Iacono et ah, 2008). Lambertian ocean surface re- 
flectance is set to 5 %. Cloud droplet effective radius is fixed 
at 15 pm, which is close to the global mean value over oceans 
observed by MODIS (King et ah, 2013). This value of effec- 
tive radius is also used to convert the MODIS visible COD to 
liquid water path used as input to RRTM-SW. Liquid cloud 
optical properties are calculated internally by RRTM. For at- 
mospheric profiles of water vapor and temperature, we use 
NCEP (National Centers for Enviromnental Prediction) R1 
reanalysis data (Kistler et ah, 2001) averaged both zonally 
and annually. Our sensitivity tests indicate that the shortwave 
DRE of ACA is largely insensitive to cloud effective radius 
or atmospheric profiles. 

3.3 Cloud optical depth correction 

As noted in previous studies (Coddington et ah, 2010; e.g., 
Haywood et ah, 2004), when a cloudy MODIS pixel is con- 
taminated by overlying light-absorbing aerosols the COD re- 
trieval is generally biased low. We have developed a fast 
COD correction scheme to account for the COD retrieval bias 
due to ACA in our DRE computation, which is illustrated 
in Fig. 2. This scheme requires both the cloud reflectance 
LUT for clouds without ACA, for which we use the MODIS 
operational LUT, and clouds with ACA, for which we use 
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the one developed by Meyer et al. (2013). In the operational 
MODIS retrieval, the reflectance LUT of cloud without ACA 
is used to interpret the reflectance of all clouds, including 
those affected by ACA. Based on this fact, we first infer the 
“observed” cloud reflectance (after atmospheric correction) 
by interpolating the reflectance LUT of cloud without ACA 
corresponding to the biased COD. Then, we use the “ob- 
served” cloud reflectance and ACA-affected LUT (derived 
based on CALIOP AOD) to determine the corrected COD. 
This COD correction process is performed for every combi- 
nation of COD bin in p(r c ) and AOD bin in /?(r a ).In the final 
step we resample the corrected CODs to obtain the corrected 
p(r c ). 

It should be noted that because different aerosol types 
may have different impacts on MODIS COD retrievals, the 
above COD correction process is aerosol-type dependent. In 
this study, we use light-absorbing aerosols as example to il- 
lustrate our method and for validation purposes we use the 
aerosol model developed by Meyer et al. (2013) for the devel- 
opment of LUTs for DRE computation and COD correction. 
However, the LUTs can be easily extended to other aerosol 
models. In fact, as part of ongoing research, we are extend- 
ing our LUTs to include all six operational CALIOP aerosol 
models as described in Omar et al. (2009). 

3.4 Uncertainty analysis 

Several recent studies suggest that the current operational 
CALIOP product tends to underestimate the above-cloud 
AOD. Meyer et al. (2013) found that the daytime CALIOP 
AOD retrievals are systematically smaller than the nighttime 
retrievals, probably due to the daytime solar background is- 
sue. In light of this finding, Meyer et al. (2013) increased the 
CALIOP AOD retrievals by a factor of 1.5 to account for the 
impact of potential AOD bias on DRE of ACA. A more re- 
cent case study by Jethva et al. (2014) suggests that CALIOP 
ACA AOD retrievals are biased low by a factor of 5 or even 
more compared with other retrievals, although the generality 
of this finding needs to be further tested with larger samples. 
While a rigorous analysis uncertainty analysis of CALIOP 
AOD product is beyond the scope of this study, it is nev- 
ertheless reasonable to assume that the current CALIOP re- 
trievals provide a lower limit to the ACA AOD. In the uncer- 
tainty analysis presented in the next section, we carry out two 
sensitivity tests to estimate the potential impacts of CALIOP 
AOD bias on DRE computation. We multiply CALIOP AOD 
values by a factor of 1.5 in the first test following Meyer et 
al. (2013) and by a factor of 5 in the second as suggested in 
Jethva et al. (2014). 

Once the magnitude of the uncertainties in the input data 
is prescribed, the consequential impact on DRE can be eas- 
ily estimated in our method as follows. First, in addition to 
the p( r a ) based on the original CALIOP data, we also de- 
rive the perturbed PDF p( r a ) by perturbing the original data 
according to predefined uncertainties (i.e., by increasing the 
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original values by a factor of 1.5 or 5). Then, the impact of 
input uncertainty on ACA DRE can be estimated by com- 
paring the DREs computed with the original vs. perturbed 
PDF (i.e., (DRE) aca vs. (DRE) aca ). Note that (DRE) ACA 
and (DRE) aca can be obtained in a single computation be- 
cause they both represent integrals over DRE(r c ,r a ), only 
with different weights. In this regard, our method is much 
more efficient than the pixel-by-pixel method, in which un- 
certainty must be estimated by perturbing individual pixels. 

4 Implementation and validation of new DRE 
estimation scheme 

Each year during austral winter, dry season biomass burn- 
ing activities throughout southern Africa inject large amounts 
of smoke into the troposphere (Eck et al., 2003; Ichoku et 
al., 2003; Myhre et al., 2003). Prevailing easterly winds dur- 
ing this season often transport the smoke to the west, off the 
continent, over the ocean, where extensive marine boundary 
layer clouds persist for most of the year. Under the descend- 
ing branch of the Hadley cell, the air mass above the bound- 
ary layer is quite dry. Due to the lack of efficient wet scav- 
enging, the transported aerosol layers can remain suspended 
in the atmosphere for days, creating a near-persistent smoke 
layer above the stratocumulus deck over the southeastern At- 
lantic Ocean (Chand et al., 2009; Devasthale and Thomas, 
2011; Keil and Haywood, 2003; Wilcox, 2010). 

To validate our method, we have compared the DRE 
of above-cloud light-absorbing aerosols in this region with 
pixel-level computations from Meyer et al. (2013). Figure 3a 
shows the seasonal mean (August/September 2007-2011) in- 
stantaneous TOA DRE of above-cloud smoke and polluted 
dust based on the pixel-level computations from Meyer et 
al. (2013). Figure 3b shows the corresponding instantaneous 
TOA aerosol radiative forcing efficiency (RFE) defined as 
the DRE per unit AOD. The DRE and RFE results com- 
puted using our method described in the previous section are 
shown in Fig. 3c and d, respectively. Evidently, both DRE 
and RFE computed using our new method agree closely with 
the pixel-level computations. Figure 4 shows the meridional 
mean DRE and RFE for the region using the results in Fig. 3. 
Not surprisingly, the outcomes of the two methods are almost 
identical. Note that the CODs used in the computations for 
Figs. 3 and 4 are directly from the MODIS products without 
COD correction. We have also compared the DRE and RFE 
from the two methods using the corrected COD and achieved 
again very good agreement (not shown because of close re- 
semblance to Fig. 3 and Fig. 4). The seasonal and regional 
mean DRE and RFE, based on the corrected COD, from 
the pixel-level computation method in Meyer et al. (2013) 
are 6.6 W m -2 and 56 Wm -2 AOD -1 , respectively (see Ta- 
ble 2). The corresponding values from our new method are 
6.4 Wm -2 and 53.8 Wm -2 AOD -1 , respectively. 
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Figure 3. (a) Seasonal mean (August/September 2007-2011) in- 
stantaneous TOA DRE of above-cloud smoke and polluted dust 
based on the pixel-level computations from Meyer et al. (2013); 
(b) seasonal mean instantaneous TOA aerosol RFE (i.e., DRE per 
AOD) from Meyer et al. (2013); (c) same as (a), but based on the 
new method; (d) same as (b), but based on the new method. Scatter 
plot of new method vs. Meyer et al. (2013) grid-mean (e) DRE and 
(f) RFE. 


Table 2. Regional and seasonal mean values of instantaneous DRE 
and RFE based on the pixel-level computation and the new method. 



DRE [Wm -2 ] 

RFE [Wm -2 


Bias adjusted 

AOD -1 ] 


(unadjusted) 

Bias adjusted 
(unadjusted) 

Pixel computation 

6.6 (5.92) 

56.0 (50.3) 

New method 

6.4 (5.77) 

53.8 (50.2) 


As previously mentioned, to estimate potential bias in 
CALIOP ACA AOD retrieval on our DRE computation, we 
carried out two sensitivity tests. We increased CALIOP AOD 
values by a factor of 1.5 in one test following Meyer et 
al. (2013) (hereafter referred to as “xl.5 test”) and by a fac- 
tor of 5 in another as suggested in Jethva et al. (2014) (“x5.0 
test” hereafter). In both cases, we corrected the MOD1S COD 
retrievals based on the scaled AOD. The regional and sea- 
sonal mean DRE of ACA increases, from the original value 
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Figure 4. Meridional mean DRE and RFE for the region based on 
the results in Fig. 3. Lines with cross symbol correspond to pixel 
computations from Meyer et al. (2013). Lines with square symbol 
correspond to results based on the new method. 

of 6.4 to 9.6 Wm -2 in the xl.5 test and to 30.9 Wm -2 in 
the x5.0 test. We have to note that this is a very rough 
estimate. Nevertheless, the DRE based on the xl.5 scal- 
ing of CALIOP AOD seems to agree reasonably with the 
value, 9.2 ± 6.6 W m -2 , reported in an independent study by 
Wilcox (Wilcox, 2012). Interestingly, the scaling of AOD 
has little impact on RFE in both cases (53.1 W m -2 AOD -1 
in the xl.5 case and 51.2 W m -2 AOD -1 in the x5.0 case), 
which apparently suggests a near-linear relationship between 
DRE and AOD as also noted in Meyer et al. (2013) and 
Wilcox (2012) (see his Fig. 5). 

In summary, as shown clearly in Figs. 3 and 4 and Table 2, 
the DRE inferred from our new method agrees very well with 
the pixel-level computations. Furthermore, the difference be- 
tween the two methods is much smaller than, for example, 
the uncertainty associated with CALIOP retrieval biases. 

It is worthwhile to clarify again that the results shown 
in Fig. 3 are seasonal mean instantaneous DRE at A-Train 
crossing time (13:30 LT) based on CALIOP above-cloud 
AOD and corrected Aqua-MODIS below-aerosol COD re- 
trievals. Moreover, the aerosol model described in Meyer et 
al. (2013) is used in this study. All these factors should be 
considered when comparing the results in this study with 
those in other studies (e.g., Chand et al., 2009; Graaf et al., 
2012; Wilcox, 2012). For example, Chand et al. (2009) used 
CALIOP in combination with Terra-MODIS observations to 
compute the DRE over the southeastern Atlantic region. It is 
known that low-clouds in this region have a strong diurnal 
cycle driven by solar cloud absorption (Bergman and Salby, 
1996; Rozendaal et al., 1995; Wood et al., 2002). As a re- 
sult the cloud properties observed by Terra MODIS can be 
significantly different from those observed by Aqua MODIS 
in this region, which could lead to different DRE even if the 
same method was used. 
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Figure 5. (a) Annual mean AOD (at 550 nm) of above-cloud light- 
absorbing aerosols (i.e., smoke and polluted dust) derived from 
4 years (2007-2010) of the CALIOP 5 km aerosol and cloud layer 
products, (b) Annual mean below-aerosol COD derived from the 
MODIS daily level-3 COD-CTP joint histogram, (c) Annual mean 
instantaneous TOA DRE of above-cloud light-absorbing aerosols 
derived using the new method. 


5 Summary and discussion 

Recent advances in satellite-based remote sensing, in par- 
ticular the launch of the space-borne lidar CALIOP, have 
provided an unprecedented opportunity for studying the 


radiative effects of ACA. However, the methodologies used 
in recent studies for computing the ACA DRE appear to 
be either oversimplified (e.g., Chand et al., 2009; Oikawa 
et al., 2013) or too cumbersome (e.g., Meyer et al., 2013). 
This paper describes a novel method recently developed for 
computing the shortwave DRE of above-cloud aerosols over 
ocean. Our method has several unique features compared 
to previous methods: (1) it takes subgrid-scale cloud and 
aerosol variation into account in DRE computations, simi- 
lar to Meyer et al. (2013); (2) it treats the overlap of aerosol 
and cloud rigorously by utilizing the joint histogram of COD 
and CTP in the MODIS level-3 cloud product; (3) it re- 
lies on grid-level cloud statistics (i.e., COD-CTP joint his- 
togram), instead of pixel-level products, and utilizes precom- 
puted look-up tables for ACA DRE computations, thus mak- 
ing it much more efficient than pixel-level computations. As 
shown in Figs. 3 and 4 and Table 2, DRE computed using 
our method agrees well with the pixel-level computations of 
Meyer et al. (2013). 

In addition to the southeastern Atlantic region, we have 
recently begun investigating the DRE of above-cloud light- 
absorbing aerosols for the global ocean. Some preliminary 
results are shown in Fig. 5. We first derived the daily grid- 
level statistics of above-cloud AOD and below-cloud COD, 
as well as the corresponding ACA DRE, using the method 
described above and then aggregated the daily means to an- 
nual mean. The temporal aggregation is weighted by the 
number of ACA pixels in each day during the period 2007- 
2010. For example, the annual mean ACA DRE in Fig. 5c 
is aggregated from the daily mean based on the following 
equation: 


E N i- (dre,} aca 

(DRE) aca = — , i= day 1,2,3..., (3) 

i 

where (DRE,) A c A is the spatially averaged instantaneous 
ACA DRE in each day averaged over ACA pixels. A, is 
the number of ACA pixels in the grid box in each day, 
and (DRE) aca is the annual mean instantaneous ACA DRE 
shown in Fig. 5c. Figure 5a shows a global map of the 
annual mean 550 nm AOD of above-cloud smoke and pol- 
luted dust derived based on 4 years (2007-2010) of CALIOP 
aerosol and cloud layer products. Similar to Devasthale and 
Thomas (2011), we note several “hotspots” of ACA over 
the southeastern Atlantic, the east-central Atlantic off the 
western coast of Saharan Africa, the Arabian Sea, and the 
North Pacific Basin off the coast of eastern Asia. It is in- 
teresting to note that the ACA AOD over the east-central 
Atlantic and Arabian Sea is noticeably larger than that over 
the southeastern Atlantic and North Pacific Basin. Figure 5b 
shows the annual mean below-aerosol COD derived from 
the MODIS daily level-3 cloud product after the correction 
of above-cloud AOD contamination using the method de- 
scribed in Sect. 3. A notable feature in the figure is that the 
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below-aerosol COD over the North Pacific Basin is signif- 
icantly larger than that over other ACA regions. Figure 5c 
shows the annual mean shortwave DRE at TOA aggregated 
from daily values due to ACA smoke and polluted dust over 
the global ocean. It is intriguing to see that the DRE of ACA 
over the North Pacific Basin is significantly larger than that 
over the southeastern Atlantic, which is in turn larger than 
the DRE over the east-central Atlantic and the Arabian Sea. 
In fact, some negative DREs are observed in the latter two 
regions. This is probably due to the COD of below-aerosol 
clouds being too thin (Fig. 5b) over these regions to have sig- 
nificant radiative effect, so that the radiative effect of ACA is 
close to that of clear skies (i.e., negative). This is interesting 
because the above-cloud AOD over these regions is actually 
larger, while the below-aerosol COD over these regions is 
smaller, compared to their counterparts over the southeast- 
ern Atlantic and North Pacific Basin. Therefore, the prelimi- 
nary results seem to suggest that the DRE variability of ACA 
is modulated by COD, rather than AOD, although it should 
be noted that we have focused only on the light-absorbing 
aerosols, i.e., smoke and polluted dust, and assumed the same 
narrowband scattering properties for them as in Meyer et 
al. (2013). Further research is needed to study the impact of 
aerosol type and scattering properties on the spatiotemporal 
variation of DRE on a global scale. Nevertheless, the prelim- 
inary results shown Fig. 5 clearly demonstrate the usefulness 
of our new method for global studies. 

It should be noted that this study, and previous ones us- 
ing CALIOP observations (e.g., Chand et al., 2008; Meyer et 
al., 2013; Oikawa et al., 2013), are limited by the capabili- 
ties of CALIOP. Arguably, some aerosols exist above every 
cloud. However, not all ACA can be detected by CALIOP 
due to its inherent limitations. Some ACAs are simply too 
optically thin to be detected, though their radiative effects 
are also expected to be small. Other situations may also be 
possible. For example, a confined aerosol layer has larger 
volume backscatter than a vertically stretched layer, even if 
the total aerosol amounts are the same, and therefore is more 
easily detected by CALIOP. Passive sensors, however, are 
less affected by the vertical distribution of ACA because they 
observe column-integrated scattering by aerosols. Recently, 
several novel techniques have been developed to detect and 
retrieve ACA properties using passive sensors. Waquet et 
al. (2009) developed a method based on multi-angular po- 
larization measurements from POLDER to retrieve the AOD 
of above-cloud smoke. This method has recently been ex- 
tended to include both smoke and dust aerosols (Waquet et 
al., 2013a). Most recently, Jethva et al. (2013) demonstrated 
the ability of a color ratio method to retrieve the above-cloud 
AOD based on MODIS multiple spectral cloud reflectance 
measurements. A review of the emerging satellite -based ob- 
servations of above-cloud aerosols can be found in Yu and 
Zhang (2013). The capabilities and limitations of the pas- 
sive techniques need to be systematically studied through in- 
tercomparisons and comparison with CALIOP observations, 
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but they may provide a complementary perspective on ACA. 
Recall that passive imagers have much larger spatial cover- 
age than CALIOP, which makes brute force calculations of 
the DRE at the pixel level computationally expensive. In this 
regard, our method satisfies the need for efficiency of ACA 
DRE computations based on passive imager retrievals. 

As a final remark, we would like to point out that the ACA 
DRE discussed in this study is still a few steps away from 
the all-sky aerosol radiative effect ((ORE) a u_ s kyl • For a given 
grid box, the (DRE) al i_ sky can be decomposed into the sum 
of clear-sky and cloudy-sky DRE: 

(DRE) al [_ sky = 

(1 - /c) ■ (DRE> clear + f c ■ /aca ■ (DRE) ACA , (4) 

where f c is the cloud fraction, (DRE) c i ear is the DRE aver- 
aged over the clear-sky portion of the grid box, /aca is the 
fraction of cloudy pixels with ACA detected by CALIOP or 
other sensors, and (DRE) A ca * s the DRE averaged over all 
ACA-containing pixels. It is important to note an implicit as- 
sumption made in Eq. (4); that is, when a distinct ACA layer 
is not detected, the DRE of ACA is zero. Different sensors 
(or different retrieval algorithms for the same sensor) may 
have different sensitivities to ACA and therefore provide dif- 
ferent estimates of /aca and (DRE) aca- For example, one 
sensor may only be able to retrieve ACA for optically thick 
clouds. This sensor would retrieve a larger (DRE) aca- hut 
a smaller /aca, in comparison with another sensor capable 
of retrieving ACA for all clouds. Therefore, when comparing 
the ACA or all-sky DRE estimated based on different instru- 
ments or algorithms, it is important to compare both the /aca 
and (DRE) aca terms in Eq. (4). 
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